Superfluidity with disorder in a quantum gas thin film 



Sebastian Krinner, 1 David Stadler, 1 Jakob Meineke, 1 Jean-Philippe Brantut/'Jj and Tilman Esslinger 1 

department of Physics, ETH Zurich, 8093 Zurich, Switzerland 
(Dated: 20121203021300) 

We investigate the properties of a strongly interacting, superfluid gas of 6 Li2 Feshbach molecules 
forming a thin film confined in a quasi two-dimensional channel with a tunable random potential, 
creating a microscopic disorder. We measure the atomic current and extract the resistance of the 
film in a two-terminal configuration, and identify a superfluid state at low disorder strength, which 
evolves into a normal, poorly conducting state for strong disorder. The transition takes place 
when the chemical potential reaches the percolation threshold of the disorder. The evolution of the 
conduction properties contrasts with the smooth behavior of the density and compressibility across 
the transition, measured in-situ at equilibrium. These features suggest the emergence of a glass-like 
phase at strong disorder. 
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Disorder has profound effects on the properties of ma- 
terials. It can turn metals pQ, or even superconductors 
and superfluids into insulators PHI]- The evolution from 
a superfluid state to an insulator depends on the details of 
the underlying material structure, which is complex and 
often only partially known. Several models have been 
devised, accounting for various aspects of the transition 
[5]. Dirty boson models, where interacting bosons are 
placed in a random potential, have attracted particular 
attention. These models predict new phases such as the 
Bose-glass, expected as a result of disorder and interac- 
tions [SHE]. 

Over the last years, ultracold gases in random poten- 
tials have emerged as an important tool to investigate 
the interplay between disorder and superfluidity. They 
offer comparatively simple systems controlled by a few 
microscopic parameters that are known and controlled 
[9]. They have been studied using observables such as 
phase coherence, response to external forces and lattice 
modulations [T0TU5| . but transport coefficients have not 
been measured, preventing a direct comparison with con- 
densed matter systems. 

In this letter, we present measurements of the DC 
resistance of a disordered, strongly interacting gas of 
6 Li 2 molecules. We use the recently demonstrated two- 
terminal setup [HI [T7] to measure the resistance of a 
quasi two-dimensional thin film of molecules as a func- 
tion of disorder strength, in a manner directly analogous 
to solid-state experiments. By comparing the resistance 
of the strongly interacting gas of molecules to that of 
a weakly interacting Fermi gas in the same potential 
landscape, we observe a breakdown of superfluid flow 
when the percolation threshold of the disorder reaches 
the chemical potential of the gas. In addition, we mea- 
sure the equilibrium density and compressibility of the 
disordered gas, which do not display large variations up 
to very strong disorder, a feature expected for a glass 
transition. 

The experimental setup is based on our previous work 
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FIG. 1: Sketch of the experimental configuration. A: top 
view. B: side view. Two 6 Li2 Feshbach molecules reservoirs 
(grey) are connected by a narrow, quasi-2D channel, created 
by a laser beam (in blue) at 532 nm with a l/e 2 -waist of 30 /im 
along y, having a nodal line at the center. When a number 
imbalance between the reservoirs is present, a molecule cur- 
rent sets in through the channel. A speckle pattern (green) 
of variable intensity is projected onto the channel along the 
2-axis through a microscope objective (not shown). 

[Trj] . and is sketched in figure [TJ We prepare a cloud 
of molecules formed during evaporative cooling of a 6 Li 
gas in a mixture of the two lowest hyperfine states at a 
magnetic field of 751 G. At this field, the s-wave scatter- 
ing length is 3545 ao, where ao is Bohr's radius, and the 
binding energy of the molecules is Ef, = 2.3 /iK [TS]. The 
elongated cloud of 1.03(5) 10 5 molecules [TO] is separated 
at its center into two reservoirs by a repulsive potential 
created with a laser beam having a nodal line at its cen- 
ter. This creates a quasi- two-dimensional channel which 
connects the two reservoirs. The trap frequency along 
the tightly confining z-axis at the center of the channel 
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is uj z — 27r • 6.1 kHz, thus huj z <C Ej,. 

Without the channel, the chemical potential of the 
molecular gas is fi = 550(70) nK, deduced from the par- 
ticle number and trap frequencies [3D] . The channel cov- 
ers a region of 30 fjjn along y, small compared to the 
size of the cloud (~ 300 /im) . Therefore we assume that 
the overall chemical potential remains constant in the 
presence of the channel. Accordingly, the local chemical 
potential in the channel is fi = /in — \hjj z ~ 400(70) nK. 
With fj, ~ lAfkjJz, the gas predominantly populates the 
lowest vibrational state. The inter-molecule scattering 
length is a m — 2100 an [2T], which yields an interaction 
parameter for the 2D scattering problem \f%Txa m jl z ~ 
1.1, with l z being the harmonic oscillator length in the 
tightly confined direction [55] . The molecular gas is thus 
in the strongly interacting regime, where the scattering 
amplitude depends on momentum |23j . 

The disordered potential is realized by projecting 
a speckle pattern onto the channel through a high- 
numerical aperture microscope. A gaussian fit to the en- 
velope of the pattern yields a 1/e 2 radius of 36 /im. The 
average disorder strength V at the center of the pattern 
is calculated from the beam envelope and the laser power. 
The correlation length a of the disorder, defined by the 
1/e radius of the autocorrelation function is 290(90) nm. 
Associated with this length scale, we introduce a corre- 
lation energy E a — =0.48 /iK, where m is the mass 
of a molecule. This scale separates the quantum regime 
V <C E a from a classical regime V > E a [H] . 

It is instructive to compare these scales to that of 
a Bose-Einstein condensate (BEC) of molecules. The 
chemical potential in the channel is independent of the 
disorder, and fixed by the unperturbed particle reser- 
voirs. The size of a molecule, given by the interatomic 
scattering length, is 190 nm, comparable but smaller than 
the correlation length of the disorder, so we expect the 
molecules to remain bound for increasing disorder. The 
healing length of the molecular BEC, associated with the 
chemical potential, is typically £ ~ 230 nm, of the order 
of the correlation length of the disorder. Thus, we ex- 
pect the disorder to affect the many-body physics at the 
microscopic level [14] . 

We investigate the conduction properties of the dis- 
ordered gas by inducing a chemical potential difference 
between the two reservoirs [THl E] • During evaporative 
cooling, a fixed relative atom number imbalance is cre- 
ated between the left and right reservoir. We observe the 
evolution of the atom number imbalance while a quasi- 
stationary current flows through the channel. We mea- 
sure the slope of the initial evolution, yielding the maxi- 
mum current driven by the chemical potential difference. 
Figure [2j\ shows the initial current as a function of V. 
For weak disorder, we observe a large current, limited by 
conservation of energy. We interpret this as a manifesta- 
tion of superfluidity, as expected for a BEC of molecules 
[25j . With increasing disorder, the current quickly de- 
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FIG. 2: Conduction properties of the disordered gas of 
molecules (red squares), and of a weakly interacting Fermi 
gas (grey dots) as a function of disorder strength normalized 
by fi = 400 nK, the chemical potential of the molecular gas in 
the channel without disorder. A: Initial current, for a fixed 
driving (see text). B: Dimensionless resistance. The dotted 
green line represents the classical percolation threshold of the 
speckle pattern. The dash-dotted line indicates the correla- 
tion energy E a . Error bars are statistical. 



creases, and the time evolution of the number imbalance 
is well described by an exponential. In this regime, we 
fit a decay time tbec = RC, where R is the resistance 
of the channel and C = dN/dfi is the V- independent 
compressibility of the reservoirs. Figure [2j3 shows r — 
RCu) yi where uj v is the underlying trap frequency of the 
reservoirs along the transport direction |17j . For weak 
disorder, we observe a fast increase of r with disorder 
strength, up to V / fi ~ 2, and a correspondingly fast de- 
crease of current. For strong disorder, the increase of 
resistance shows a weaker but again exponential depen- 
dency on the disorder strength. The highest resistance 
plotted represents a slow down of the transport by about 
three orders of magnitude compared to the case without 
disorder. It represents the lowest current and highest 
resistance we can measure. 

To disentangle the effects of superfluidity from the sin- 
gle particle effects caused by disorder, we repeat the ex- 
periment in the same trap and disorder configuration 
with a weakly interacting Fermi gas (WIF) of 9.3(4) 10 4 
atoms per spin state, prepared at a magnetic field of 
475 G, where the scattering length is — 100 ao- The cur- 
rent and resistance are shown as grey dots on figure [2] 
They show a single exponential dependency and a finite 
resistance at low disorder strength, corresponding to the 



3 



contact resistance of the channel [TB] . This behavior con- 
firms that the disorder strongly affects transport already 
at the single particle level [26H28] . 

We then extract a parameter B = Rbec/RwiFi which 
measures the effect of pairing and superfluidity com- 
pared to the ideal Fermi gas case, keeping the geom- 
etry and potential landscape unchanged, ft accounts, 
at least qualitatively, for the variations of resistance in 
a non superfluid case, without relying on a model for 
the complicated single particle dynamics. This proce- 
dure resembles the use of the normal state resistance 
as a reference in the physics of superconductors. We 
evaluate B — -^tbec/twif where K is the ratio of the 
compressibility of the reservoirs |29j . which does not 
depend on disorder strength. Modeling the reservoirs 
as harmonically trapped zero temperature Fermi gases 
and weakly interacting Bose-Einstein condensate, we es- 
timate K ~ 0.57. 
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FIG. 3: B parameter (see text) as a function of disorder 
strength. The green dotted line represents the percolation 
threshold of the speckle potential, and coincides with the 
change of behavior from superfluid to single-particle trans- 
port. The dash-dotted line shows the correlation energy E a . 
Error bars represent statistical errors. 

Figure [3] shows the evolution of B with disorder 
strength, obtained by the ratio of the decay times for the 
molecules and the linearly interpolated data of the WIF. 
For strong disorder, it remains constant. The exact value 
should depend on the chemical potentials of the reservoirs 
for the two cases, but it should reflect the zero temper- 
ature behavior of conductance [3D]. This indicates that 
transport in the strongly disordered regime is dominated 
by single atom or molecule physics (26j . In contrast, for 
weak disorder, B varies very quickly, incompatible with 
simple diffusion 26J . B displays a sharp change of behav- 
ior for intermediate disorder strength, at the point where 
the chemical potential of the reservoirs reaches the perco- 
lation threshold of the disorder V t — 1.92// [3TJ|32] (green 
dotted line). This corresponds to the intuitive picture of 
a transition from isolated superfluid pockets to a con- 
nected superfluid allowing for superflow. It agrees with 
the concept of a Bose-glass composed of an incoherent 
ensemble of localized condensates [24;, 33-36 . However, 



contrary to the superfluid to Bose-glass quantum phase 
transition, the behavior in the strongly disordered case 
resembles that of a metal with finite conductance. Since 
the Bose-glass is smoothly connected at finite tempera- 
ture to the normal Bose gas, we propose that the finite 
resistance is a finite temperature effect. It is remarkable 
that the classical picture of percolation gives a quanti- 
tative description of our observations even though in- 
teractions are strong and tunneling through individual 
disorder grains is significant, as V t ~ 1.6E a |24) . 

We now investigate how the density and compressibil- 
ity in the channel depend on the disorder strength. We 
prepare a cloud with equal population in the two reser- 
voirs. The disorder is switched on and an absorption pic- 
ture is taken in-situ after 150 ms, allowing the channel to 
thermalize with the reservoirs [37] . We typically average 
20 of those pictures to reduce noise. We first consider the 
density observed at the center of the channel, averaged 
over a region of 18 /zm (7 /Ltm) along the y (x) direction. 
Over this region, the trapping frequency along z varies 
by less than 10 %, and the local chemical potential of the 
cloud by about 30%. The area of this region is much 
larger than the correlation length of the disorder, there- 
fore we expect the spatial average to reflect the average 
over disorder realizations. 

Figure [4]\_ presents the number of atoms per corre- 
lation area of the disorder h — rnia 1 as a function of 
disorder strength. The density smoothly decreases with 
increasing disorder up to the highest disorder strength. 
This is expected since the disorder is repulsive and re- 
duces the available phase space. The density changes by 
a factor of 5, which is small compared to the three orders 
of magnitude change in the resistance. Therefore, the 
fast variation in resistance can not be solely attributed 
to a change in density. 

We now make use of the in-situ density measure- 
ments to evaluate the compressibility. Because the trap- 
ping potential along x varies much slower than the 
correlation length of the disorder, we use the local 
density approximation to obtain the variation of the 
disorder-averaged density with chemical potential, i.e. 
the disorder-averaged compressibility dn/dy,, where n 
is the disorder-averaged density. At the center of the 
cloud, the curvature of the density distribution reads [33] 
gjS 1 = — rnw^^j, where we include in uj x the small effect 
of the envelope of the speckle pattern. We use a quadratic 
fit to the central part of the cloud to extract the curva- 
ture of the density profile, shown on figure [4p and D for 
two different disorder strength. 



Figure 



presents the reduced compressibility ^^j- 



For very weak disorder, the compressibility remains con- 
stant, equal to that of the system without disorder. For 
stronger disorder, it decreases continuously with disorder 
strength. The overall variation remains within a factor 
of ~ 5. The compressibility remains finite up to very 
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FIG. 4: Equilibrium properties of the disordered gas as a 
function of disorder strength. The green dotted line indi- 
cates the percolation threshold. The solid black line indicates 
the prediction of the Thomas-Fermi approximation. A: av- 
erage density at the center of the channel. The dashed line 
indicates the expected decrease of density from the compress- 
ibility measurement without speckle (see text). B: Reduced 
compressibility (see text) estimated from the shape of the 
cloud. Insets: density profiles in the channel (open circles) 
and quadratic fit to the central part (continuous line) used to 
extract the compressibility, for V/fi = and 1.6 (C and D). 

strong disorder where transport is frozen, suggesting a 
glassy phase. 

The zeroth order effect of a weak repulsive disorder is a 
shift of the average chemical potential equal to the mean 
disorder strength. We can use this and the measured 
compressibility to estimate the decay of the density due 
to the disorder (dashed line in figure [4j\) . The agreement 
gives us confidence in the compressibility measurement. 

We can further compare these observations to the pre- 
dictions of a zero-temperature Thomas-Fermi approxi- 
mation for a random potential |36] , shown as a black line 
in figure |4]A. (B). The free parameters are the chemical 
potential and compressibility in absence of disorder. We 
find good agreement at low disorder, where the disor- 
der mainly gives an overall shift of V of the chemical 
potential. For higher disorder, the observed density lies 
below the predictions. This is expected as h < 1 and in- 
teractions are strong, so that the fragments of BEC are 
strongly depleted. Surprisingly, the shape of the com- 
pressibility data is well reproduced. 

Our measurements have shown the difference between 
equilibrium and transport properties of a superfluid 
bosonic thin film with disorder. By tuning the inter- 



action using the Feshbach resonance, our setup could be 
used to investigate another, Fermionic scenario, where 
pairing results from a many-body process and can be 
strongly affected by the disorder SO] • 
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